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Abstract 

I investigate the time step dependence of Monte Carlo simulations 
for coordinate-spaces consisting of several patches. It is shown that a 
naive kinetic term does not necessarily converge to the same spectrum 
as a Hamiltonian calculation. Then an improved kinetic term is pre- 
sented which allows one to connect the Monte Carlo and Rayleigh-Ritz 
results of intermediate volume SU(2) gauge theory. 

1 Introduction 

There have been two approaches in the intermediate volume calculation of 
SU(2) gauge theory: The Hamiltonian approach [||], which is a Rayleigh- 
Ritz analysis of an effective Hamiltonian. and the Lagrangian approach |^ 
which is a Monte Carlo simulation of an effective Lagrangian. For a review 
of most of the work done in this area see 0]. Although both approaches 
are formally equivalent 0, the numerical data obtained does not agree. A 
study of the Lagrangian approach using smaller time steps did not improve 
the results 0. 

In both approaches the coordinate space consists of 8 compact patches, 
which correspond to equivalent vacua. Basically these patches are spheres 
and tunneling from one sphere into another one can take place via any point 
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of the boundary of the sphere. This means, that in flat space, we cannot 
form a single large patch by joining all the patches together. The idea, 
that this kind of geometry is responsible for the difference of the Lagrangian 
and Hamiltonian results, is very tempting. Indeed, the disagreement in the 
energy spectrum is still present, when one switches off the potential energy. 
One can therefore study the simpler case of a free particle moving in the 
same geometry, which allows one to use analytical methods. This is carried 
out in section 2. After having understood the reason for the disagreement, 
one can change the kinetic term in the Lagrangian approach to match the 
Hamiltonian approach. A procedure to do so and some Monte Carlo results 
are given in section 3. Finally section 4 contains some conclusions. 

2 A toy model 

The model without the potential energy factorizes into 3 wave functions. 
For the analysis attempted here, it is sufficient to describe one of these wave 
functions: Our toy model consists of a free particle moving in two spheres, 
each of radius tt. The obvious choice of coordinates are spherical coordinates 
and an integer A; = 0, 1, giving the sphere in which the particle is: 

i^^^{r,e,ip,k). (1) 

The spheres touch each other at their boundary. (One can picture this a 
two spheres occupying the same space.) For the particle to move from the 
interior of one sphere to the other one, it has to move to any point on the 
boundary r — tt where it changes to the other sphere and then moves on 
inward. 

Due to the symmetry of the problem the eigenfunctions will be either 
symmetric or anti-symmetric under interchange of the spheres. One can 
therefore describe the wave function by its wave function ip in one sphere 
and its symmetry type. 

For the anti-symmetric wave functions the boundary condition at r = tt 
is unambiguous: 

^l^ir,9,ip) = 0. (2) 

For the symmetric wave function the choice is not as clear cut. Due to the 
measure one should however look at r'0(r, 9, ip) instead of ■0(r, 9, (p). One can 
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select a basis B{b) by requiring 

d 



(3) 



for a given b. In the original problem the angular variables could not be 
observed, and the obvious choice was 6 = 0. Here, however, one is free to 
choose any b. The eigenfunctions are given by 

ipiim{r, 0, if) = jiiair)Yim{0, if) (4) 

where is the i'th zero or @. That these functions form a basis for 
any 6, follows from the orthogonality properties of the Bessel functions. 

In the Lagrangian approach, or to be more precise the Monte Carlo cal- 
culation with finite time step e P], the configurations are described by a 
spatial coordinates r{t) and a coordinate k{t) G {0, 1} which describes in 
which sphere the particle is. The measure is flat in r and one possible kinetic 
term is given by 

TUm,r{t+e),Ak) = ^(r(t)-r(t+6))^+Afc ^^' ~ '^^'^^2,2 '^^^ ^ > 

(5) 

where Afc = \k{t) — k{t + e)|. In |^ it was proposed that, for AA; = 1 one 
should calculate the geodesic distance dg between c(t) and c(t + e) and use 

= (6) 

It was argued that both kinetic terms should give the same continuum limit, 
but (ID should have smaller finite e effects. Unfortunately the geodesic dis- 
tance is the root of a quartic equation which makes analytic calculations hard. 
I shall therefore restrict a transfer matrix analysis to the Lagrangian (|^). 
Later on I shall do a numerical comparison with (P). 



2.1 Transfer matrix calculations 

One can calculate the energy spectrum for finite time step e by means of the 
transfer matrix approach. To this end let me calculate the transfer matrix 
between two wave functions ip and 

'^i>,4>= H dfldf*2i'{fl,ki)(f)*{f*2,k2)K{fl,r2,\ki- k2\), (7) 

fcl,fc2=0 



3 



where 



V2 



■exp [-eTfei„(fi,f2, Afc)] . 



(8) 



7re 



Due to the symmetry we can spht up the transfer matrix into and interior 
(To) and exterior (Ti) part: 



T(v;,0) 



{ip , (f)) + {ip , (p) , -0 and symmetric 
T'^ {ip , (p) — {ip , (p) , -0 and antisymmetric (9) 
, otherwise 



with 



T'=(^,0) = / dr-ldr-2ij{r-l)4>*{r-2)K{r-l,r2,k), 



(10) 



where tp and are wave functions of the form (^. Due to the spherical 
symmetry of the problem, wave functions with different angular momentum 
do not mix and one has 



T {ipilmi 4'jpq) — '^ij.l^lp^rnq- 



(11) 



2.2 Expansion for small time steps 

The matrix elements can be calculated as a series in e^/^. I shall restrict 
myself here to the more interesting case of the symmetric sector. Also I shall 
give results for / = only and drop the / and m references and angular 
variables from the notation. One obtains for the diagonal elements: 



1 + 



127r(7r62 -b + Tiaf) 



[ate - 6a^e - 12nb + 6) 



+ 0(e 



3/2^ 



(12) 

(Note that can be arbitrary large, that is afe is not necessary small.) For 
the off-diagonal elements, of order 1 and aj arbitrary, one has 



+ 0{e^'^), (13) 



with 



(14) 



4 



For ttj of order 1 one will need the next order in e^/^ as well, and the 
matrix element for both, and aj of order 1 is given by: 

Let us assume we can do perturbation theory around \Ef = ipii^r). The 
first order corrections to the groundstate wave function are 

(^.^) = (16) 

For small aj one has 

i^j^H) = (27r - 6)0(1) + 0(e^/2). (17) 
and for aj > 0(1 /y/e) one has 

iS/y - 6 + ™g a2(eH'-l) 

Note that for b = the corrections (0) to the wave function vanish 
exponentially, once > e, whereas for all other choices for b they only 
vanish as a~^. Thus, for any finite e we expect the correct wave function 
to obey boundary conditions (|^) with b = but the region where the 

slope of the wave function becomes compatible with this b.c. becomes more 
narrow as e decreases. In this simple argument I ignored the fact that in the 
basis b = —-^ the first 0{1/^) trial functions contribute substantially to 
the true groundstate wave function. Later on I shall give a more rigorous 
derivation of the boundary condition. 

From eq. (|1^) one sees, that the only basis were one is allowed to do 
perturbation theory (in e^^^) is 6 = 



For this choice of b the matrix element (18) reads: 



/47raf - 1 (^ea"^' - Ij 
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For ttj < 0(e ^/^) the dominant contribution comes from the 0(e^/^) term 
in (0): 

37r3/2(a2 - af)^iA7i^a^ - l)(47r2a| - 1) 

Let me calculate the sum over all matrix element squared to see whether 
perturbation theory can be apphed to this basis. From eq. (||) one has aj — 
dj-i ~ 1, and because we are only interested in the order of the correction, 
we can replace the sum by an integral. Eq. ([19D gives a contribution 

da '\ J ' =0{e'l\ (21) 
^ "'0 a2(e2" - 1) 



whereas eq. (|^) contributes 



^ / c^a- = 0(e). (22) 



Thus, corrections to the groundstate wave function vanish as e ^ 0, and one 
can indeed do perturbation theory about e = 0. 

One is now in a position to determine the energy spectrum Ei = — log(Tjj)/e 
of the Lagrangian. I am using the basis b = ^ and restrict myself to the 
low- lying states = 0(1). 

In addition to the diagonal element Ta one needs the second order cor- 
rection 

AV. = Y.m, = j:^^^.- (23) 

The only term contributing to 0(e'^/^) is 

4^2 2 f(2 + a2e)e-H^ -2^ 

Ay..- Z5_! L (OA) 

"""^^^ - vr2(4vr2af - 1) a|e2(l - 6"^^) ^ ^ 

For small j the contribution Al^j from (p^ vanish as and can be neglected. 
Therefore one can replace the sum over j by and integral over a to obtain 

4a2e3/2 ^oo ((2 + a2)e-5«' -2)' 
7r2 47r2a2 - 1 io 0^(1 - e"^"') ^ ^ 
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Together with the diagonal element Tu up to 0(e^/^) 
one obtains for the energy 

Eq. (D) gives Oq = 0.3710096482 . . ., and one gets for the groundstate energy: 

Eq = 0.06882408 - 0.0085092 v^. (28) 

Using the same approach in the ant i- symmetric sector one obtains for the 
energy of the lowest state 

K = \ + ^e. (29) 

Let me come back to the boundary conditions in the symmetric sector. 
We can confirm our previous analysis about the derivative at the boundary 
by calculating it in the basis b = To this end let me write 

^ log(r«(r. M)) = ■im /<">-^<7;^> . (30) 

or «-*o MVt$(7r) 

where $(r) = r\E'j(r, ^, 0). Let me restrict myself to the case / = 0. Except 
for ipQ the only relevant contribution come from large j, and one has: 



$(7r) - <l>(7r - MVe) ^ 



3/2,/47r2a2 - 1 



WK cos(au) — 1 

-— / da 



ia2 
62 



(31) 



The integral in eq. (^) is of order and its value can be ignored. One 
obtains for derivative 

$(7r) - $(7r -UyTe) _ l 

u^e^TT) ~ 2n^^^'''- ^^^^ 

Thus we arrive the previous result. Because the first correction in u is 0{u) 
(as opposed to 0{ue^^'^)), one has <l>'(r)/$(r) ^ for |r — 7r| < ^/e. 
But further away from the boundary and small enough e the wave function 
becomes compatible with the lowest function of the b = basis. 
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Table 1: Rayleigh-Ritz results. Shown are the number of basis states rib 
used in the Rayleigh-Ritz calculation, the groundstate energy Eq, the lowest 
state in the anti-symmetric sector Ea and the lowest / = 1 state Ei . All Ei 
estimates were computed with 50 basis states. 



e 


rib 


^o(e) 


Ea{e) - 


1 

2 






2.0 


10" 


-1 


50 


0.066552 


1.27270 


■ 10- 


-2 


0.308909 


1.0 


10" 


-1 


50 


0.066858 


5.92709 


• 10- 


-3 


0.306617 


5.0 


10- 


-2 


200 


0.067272 


2.82432 


■ 10- 


-3 


0.305885 


2.0 


10- 


-2 


300 


0.067757 


1.08418 


■ 10- 


-3 




5.0 


10- 


-3 


300 


0.068256 


2.61903 


■ 10- 


-4 




2.0 


10- 


-3 


400 


0.068457 


1.03473 


■ 10- 


-4 




2.0 


10" 


-4 


600 


0.068706 


1.01993 


■ 10- 


-5 




2.0 


10" 


-5 


800 


0.068786 


1.01536 


■ 10- 


-5 





2.3 Rayleigh-Ritz and Monte Carlo results 

Instead of expanding the matrix elements Tij in a series in e^/^, they can be 
calculated numerically (ignoring terms of order exp[— 27r^/e]) and used for a 
Rayleigh-Ritz analysis. The internal elements 7^° can be integrated analyt- 
ically, while the external elements T-^ can be reduced to a one dimensional 
integral which has to be performed numerically. Due to the exponential fall- 
off of the matrix elements Tij for large j, the basis b = is best suited 
for this calculation. Depending on e, between 50 and 800 basis states were 
used. For the three largest e, the transfer matrix in the / = 1 sector was 
calculated by means of a two dimensional numerical integral. The results for 
the vacuum Eq, the lowest anti-symmetric state Ea, and the lowest / = 1 
state El are given in table |l[ From the results one can estimate the next 
term in the expansion of -E'o(e) in term of e: 

Eo{e) ^ 0.06882408 - 0.0085092v^ + 0.007e (33) 

In figure ^ the Rayleigh-Ritz results as well as the two series are shown. 

Using the data from table |l| one can calculate the mass gaps for the first 
antisymmetric state and the first / = 1 state mi. I have also calculated 
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1e-05 I ■ ' — ^ ' ■ — ^ ■ ■ — ■ ' — ^ ■ ■ — ^ 

1e-05 0.0001 0.001 0.01 0.1 1 

e 

Figure 1: Comparison of the Rayleigh-Ritz calculation (diamonds), expan- 
sion up to e^/^ (solid line) and estimated expansion up to 0(e) (dashed line) 
for the groundstate energy AEq = Eo{e) — Eo{0). 

the mass gaps for these sectors by means of Monte Carlo simulation, using 
both kinetic terms, (||) and (||). Table ^ shows the results. 

The data for the asymmetric state is also plotted in figure 0. One can 
see good agreement between the Rayleigh-Ritz calculation and the Monte 
Carlo simulation with kinetic term (^. From the data it seems however very 
unlikely that the results for the different kinetic terms converge at e = 0. In 
addition, a look at the groundstate wave function (which can be obtained 
from histograms in the MC calculation) shows a different behaviour for the 
two kinetic terms as well. This results comes as no surprise if one looks at 
a very simple case: In it was shown that for f(t) and r(t + e) equal but 
in different patches, the kinetic term (||) becomes Tkin = (25)^(1 — (5/(27r))^, 
where 5 = vr — r{t). The geodesic distance for this configuration is exactly 
(2(5)^. If one follows terms of the form through the small-e expansion, one 
finds they contribute 0(e) to the external matrix element. But this exactly 
the order which gives the zeroth order energy. 

As mentioned earlier the geodesic distance is the solution of a quartic 
equation. Even for small e, configurations which have multiple zeroes of the 
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Table 2: Monte Carlo results for the toy model. The table shows Rayleigh- 
Ritz and Monte Carlo results for rua and mi. RR denotes the Rayleigh- 
Ritz calculation, MC the Monte carlo results obtained by using the kinetic 
term and gd are the masses for the simulation with kinetic term (j^). 



e 


rria 


mi 




RR 


MC 


gd 


RR 


MC 


gd 


0.05 


0.435553 


0.436(5) 


0.336(4) 


0.23862 


0.2378(15) 


0.2672(10) 


0.10 


0.439069 


0.440(4) 


0.338(2) 


0.23976 


0.2394(10) 


0.2696(10) 


0.20 


0.446175 


0.446(2) 


0.341(2) 


0.24236 


0.2420(4) 


0.2668(6) 



0.05 0.1 0.15 0.2 0.25 

Figure 2: Results for m^. The solid line is the Rayleigh-Ritz calculation, the 
diamonds are the Monte Carlo simulation using the kinetic term (^) and the 
crosses are the results for the kinetic term (P). 
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equation are not suppressed. The influence of tliese configurations is beyond 
the scope of this article, but I would suspect that they also contribute to the 
observed energy difference. 

3 An improved kinetic term 

Instead of using a kinetic term like it is possible to use a kinetic term 
which will exactly reproduce the Hamiltonian result for the toy model. Again 
one would be free to tailor the kinetic term to any boundary condition (^. In 
order to make connection to the Rayleigh-Ritz results of 0, I choose b = 0. 
1 make the Ansatz: 

„ 1 oo 

s=0 j,l=0 

(34) 

where c is the angle between ri and r2 and aoij and auj are chosen according 
to @ and (H) respectively. For large enough e one can stop the sum at 
some values Imax and jmax without introducing a large error (compared to 
the statistical error of the MC-calculation). For e down to 0.15, I found 
Imax = 100 and jmax = 40 to bc sufficient. Of course it is impracticable to 
calculate the sum at every update step of the Monte Carlo procedure. But 
one can tabulate the values and look them up during the actual calculation. 
To this end we notice that (0) depends on three continuous coordinates ri,r2 
and c. The influence from the boundary will be large only for fi and r2 close 
the boundary. Otherwise 

- eLkinin, f2, 0) ^ -^('^i + ^2 - 2rir2c), (35) 

and configurations with Ak = 1 and r*! or r2 away from the boundary are 
strongly suppressed. Therefore I chose to tabulate L^in as a cubic spline 
interpolation S{ri,r2,c) in 3 dimensions, which is capable of describing (^) 
exactly. For spline interpolation the gridpoints do not have to be equally 
spaced (for a discussion of spline interpolation see for example 0). A rough 
estimate suggest that the influence of the boundary, i.e. the deviation from 
the quadratic form ( pS]) is of the order of exp[— (tt — r)^/(2e)]. Therefore 
the gridpoints (in both r directions) are chosen according to the distribution 
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Table 3: Monte Carlo results using the improved kinetic term. 



state 


exact 


MC 


Ai 


1.0 


0.9985(15) 


A2 


3.0 


2.998(4) 


e 


0.375 


0.373(2) 


P 


0.25637.. 


0.2563(6) 



function 

F(r) =pr + erf(7r/V2i) - erf((7r - r)/V2i), (36) 

where j3 is chosen to minimize the maximum errorQ. For the spline interpo- 
lation one needs to store for each gridpoint eight values: 

^ ^ d^S d^S d^S d^S 

' drf ' ' dc^ ' drldrl ' drfdc^ ' drldc^ ' drldrldc^ 

The eight quantities (|37D can be calculated from (|3^ . (The derivatives 
at the boundary which are needed to calculate (pT]) can also be obtained 
from eq. (|3^) ). Once the spline table is set up the calculation of the kinetic 
energy takes about 200 floating point operations. 

For the simulations I used a 50'^ grid. Due to the symmetry under ri <-> r2 
the size of the table is 8 x 2 x 50'^/2 numbers, which translates to 8 Mbytes 
when using IEEE double precision. The placement of the r-gridpoints are 
shown in figure ^ In figures ^ and |^ an estimate of the error | exp[— eS*] — 
exp[— eL]| of the spline approximation to the Boltzmann weight is shown. 

For all e used in the simulation the maximum absolute error was about 
10^^ and can be neglected against the statistical error of the simulation. Due 
to the Ansatz ([3^ ) the Monte Carlo calculation of the free particle should 
reproduce the Hamiltonian results for any value of e. In table (^ data for 
e = 0.6 is given. Smaller e give the same results within errors, but need more 
CPU-time to obtain the same accuracy. 

However, when one returns to the original problem there will be O(e^) 
corrections, due to the potential energy term, and one has to estimate the 

^if j3 is too small the interpolation starts to oscillate away from the boundary 
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i ' ' ' ' ' ' ' ' ' 1 

5 10 15 20 25 30 35 40 45 50 

support point 

Figure 3: Location of support points for spline interpolation. Shown are the 
locations for e = 0.6 (diamonds) and e = 0.15 (crosses). 



1.5e-08 - 




Figure 4: Absolute error of sphne interpolation to the Boltzmann weight for 
both coordinates in the same sphere 
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Figure 5: Absolute error of spline interpolation to the Boltzmann weight for 
the coordinates in different spheres 



finite time step errors. It is therefore necessary to simulate the Lagrangian 
at different time steps e. To clarify the notation used for the results let me 
write the Lagrangian in its original form: 

~eL = \{g-^ + ai)^ + ~eV{r). (38) 



One can connect this form to a free particle and eq. (|3^) by writing the 
kinetic term as 

l(Af)2 ^ ^ 

ITk^n = 2^^- (39) 

This gives for the effective time step 

e = a + ai j , (40) 

which has to be used to calculate the mass gap from the correlation function. 
The kinetic energy term, however, has to be computed by using a for the time 
step in eq. (|3^) . Monte Carlo simulations were performed for a = 0.15, 0.3, 0.6 
and g = 1.4,1.8,2.0. The results for various states are given in table (H), 
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which also contains an extrapolation down to a = 0. One finds excellent 
agreement for most states, indicating that the discrepancy of the previously 
published results is entirely due to the kinetic term. 

4 Conclusions 

The treatment of non-trivially connected coordinate spaces using Monte 
Carlo methods needs some care. There are basically two effects which hinder 
the MC approach with a naive kinetic term. Firstly one encounters correc- 
tions of the order of the square root of the time step, which require very 
small time steps and therefore large amounts of CPU time. Secondly, and 
more severe the conditions at the coordinate patch connection are not under 
control. Usually one would argue that the Monte Carlo approach will choose 
the correct boundary conditions. But if the connection of the patches is not 
simple, this argument is on very weak footing, and the boundary condition 
(and spectrum) preferred by the simulation depends on details of the kinetic 
term. The problem can be overcome by constructing a kinetic term which 
will match predetermined conditions at the boundary. In the case studied 
here, the term could be tabulated and used efficiently in a MC calculation. 
The improved kinetic term has no finite time step correction, and one is left 
with O(e^) correction from the potential energy. 

Using the improved kinetic term, I was able to reproduce the Rayleigh- 
Ritz results of |^ with the Lagrangian approach. Note that both approaches 
have the same degree of freedom in choosing the boundary conditions. In 
the original Lagrangian approach the freedom of choice is not as obvious, 
because it is hidden non-trivially in the choice of sub-leading orders of the 
kinetic term. 

Let me remark at this point that the transfer matrix approach of is 
unrelated to the approach presented here, that is van Baal's result of small 
finite time step corrections does not contradict my findings. 
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Table 4: Masses for SU(2) small volume calculation. The values for a = 
0.15,0.3,0.6 are measurements. The a = column is the extrapolation to 
a = assuming an dependence of the mass. Cont. is the result from 
the Hamiltonian approach For some states I was unable to extract the 
mass for a = 0.15, the corresponding entries are left blank, e, €2 and are 
states with 1, 2 and 3 units of electric flux. Af , ,T2' , Ai are IR of the 
cubic group, is a A^ IR with 3 units of electric flux and is a B2 
representation of with one unit of electric flux. 



g = lA 


state 


a=0.6 


a = 0.3 


a = 0.15 


a = 


cont. 


e 


0.2270(3) 


0.2034(6) 


0.1973(10) 


0.1954(6) 


0.195 




0.480(3) 


0.430(5) 


0.418(3) 


0.414(2) 


0.416 


es 


0.787(6) 


0.709(9) 


0.689(10) 


0.683(8) 


0.686 


At 


1.35(6) 


1.40(5) 


1.52(3) 


1.51(3) 


1.509 


E+ 


1.323(6) 


1.337(10) 


1.363(10) 


1.356(8) 


1.354 




2.552(15) 


2.52(3) 




2.51(4) 


2.523 




0.859(3) 


0.950(10) 


0.999(7) 


1.002(7) 


1.007 




5.42(15) 


5.5(2) 


5.52(5) 


5.53(5) 


5.479 




6.08(10) 


6.1(2) 




6.1(3) 




(7 = 1.8 


state 


a=0.6 


a = 0.3 


a = 0.15 


a = 


cont. 


e 


0.523(1) 


0.507(2) 


0.506(3) 


0.503(2) 


0.503 




1.133(4) 


1.103(6) 


1.098(15) 


1.094(7) 


1.095 


es 


1.92(2) 


1.86(3) 


1.86(3) 


1.85(3) 


1.859 


At 


2.43(5) 


2.50(5) 


2.56(10) 


2.54(6) 


2.509 


E+ 


2.27(1) 


2.254(10) 


2.27(2) 


2.255(10) 


2.258 




4.46(4) 


4.39(5) 


4.35(7) 


4.36(5) 


4.350 




1.338(5) 


1.380(5) 


1.388(10) 


1.393(6) 


1.386 


A-, 


5.92(4) 


5.98(3) 


6.02(4) 


6.01(3) 


5.987 




8.6(2) 


8.4(1) 


8.3(2) 


8.3(2) 




(7 = 2.0 


state 


a=0.6 


a = 0.3 


a = 0.15 


a = 


cont. 


e 


0.7442(10) 


0.7312(15) 


0.727(3) 


0.727(3) 


0.728 


£2 


1.621(7) 


1.588(10) 


1.577(10) 


1.575(8) 


1.586 


es 


2.72(2) 


2.69(3) 


2.67(3) 


2.67(3) 


2.672 


At 


3.15(15) 


3.30(10) 


3.1(2) 


3.27(15) 


3.177 


E+ 


2.879(10) 


2.884(15) 1' 


5 2.84(3) 


2.87(2) 


2.864 




5.54(10) 


5.56(5) 


5.49(7) 


5.53(5) 


5.513 


ht 


1.642(5) 


1.657(10) 


1.62(3) 


1.656(15) 


1.660 


Ar 


6.05(2) 


6.08(3) 




6.09(4) 


6.113 


at 


10.2(3) 


10.3(2) 


10.1(2) 


10.2(2) 


9.79 
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